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Abstract: The present contribution deals with a multi-compartment model 
simulating nymphal growth of hemimetabolous insects. In general, a horwmonal growth 
rate (Wn) can be distinguished from an intrinsic growth rate (Wi), whereby the second 
factor includes the increase of nymphal size due to food consumption. Besides this basic 
assumption the theoretical approach also considers the possible influence of external 
factors (environmental temperature, available food resources, circadian cycle, 
intraspecific competition) on insect growth. Model application and validation were 
carried out for Teleogryllus commodus WALKER 1869, where also experimental data of 
nymphal development are available. Preliminary modeling results correspond very well 
with comparable experimental data (r: 0.97-99). Basically, a significant dependence of 
nymphal growth on environmental temperature can be recognized. 


Key words: Nymphal growth, hemimetabolous insect, multi-compartment 
model, growth rate, cricket, Teleogryllus commodus. 


Introduction 


Since the legendary work of J. B. S. HALDANE (1927) research of growth physiology 
within the class of insects has experienced a considerable accretion. Until the middle of 
the last century, insect science favored the hypothesis that growth and maturation of a 
given species can be exclusively regarded as the consequence of food intake (BAKKER 
1959, HOFFMANN 1985). Currently, it is considered a proven fact that holometabolous 
insects are able to temporarily disconnect their larval growth from direct food intake, 
because they can hark back to respective energy reserves stored in the body. This process 
is mainly controlled by insect hormones (JH, EH, PTTH) and the cellular signal cascades 
induced by them (CONLON & RAFF 1999, STERN 2003). In the case of the fruit fly 
Drosophila melanogaster organ and body growth is additionally influenced by so-called 
DILPs (Drosophila insulin-like peptides), whereby quantitative production of these 
hormones shows a remarkable dependence on the amount of ingested food and occurs in 
small cell agglomerations of the insect brain. Intensity of larval growth correlates with 
the concentration of DILPs in the hemolymph. According to a recent theory these 
hormones regulate the absorption of nutrients from the hemolymph into the the cells of 
the growing structures (BAKKER 1959, NIJHOUT et al. 2006, SHINGLETON et al. 2007). 


With regard to hemimetabolous insects higher numbers of larval (nymphal) instars 
associated with greater time spans of nymphal development can be observed. Most 
cricket species are characterized by a duration of larvogenesis ranging from 50 to 100 
days, whereby the organism strides through 8 to 14 instars (Fig. 1; WEBER & WEIDNER 
1974, CHAPMAN 1998, SHINGLETON et al. 2008, STURM 1999, 2002a, 2003a, 2016a). As 
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found by comprehensive experimental studies, each nymphal instar of a hemimetabolous 
insect is distinguished by individual growth rates being associated with particular 
concentrations of various growth hormones (MERKEL 1977, BEHRENS et al. 1983, 
HOFFMANN 1985). 
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Fig. 1: Development of hemimetabolous insects (e.g., crickets) with a predefined number of 
nymphal instars. Duration of larvogenesis among other depends on external factors such as 
environmental temperature. 


Theoretical models dealing with the larval or nymphal growth of insects are chiefly 
founded upon specific cycles of food intake and hormone production. It is hypothesized 
that the quality and quantity of the ingested food has a remarkable influence on the DILP 
level, establishing a relationship between physiological processes on the one hand and 
environmental conditions on the other. Previous studies could demonstrate that 
environmental temperature may be regarded as primary control factor concerning insect 
development (DAMOS & SAVOPOULOU-SOULTANI 2008, 2012). In early growth models 
this phenomenon was considered by definition of the law of total effective temperatures 
(HILBERTAND & LOGAN 1983, WAGNER et al. 1984), whereas modern approaches are 
based on the assumption that insects only perform an optimal development within a 
specific temperature interval. Outside this thermal range developmental rates are subject 
to a continuous decrease and final cessation (DAMOS & SAVOPOULOU-SOULTANI 2008, 
2012, STURM 2016a). 
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In general, linear models describing a linear relationship between temperature and insect 
development can be distinguished from non-linear approaches, where the thermal effect 
on nymphal growth is commonly expressed by polynomials of higher degree. The second 
model type is characterized by the circumstance that developmental rates become maxi- 
mized at optimum temperatures and rapidly approach the zero-value outside this prefer- 
ential thermal frame. Most current models offer a simulation of insect development in 
biophysical and biochemical terms, whereby enzymatic reactions play a superior role. 
According to the equations formulated by VAN’T HOFF (1901), ARRHENIUS (1889), and 
EYRING (VAN DE HAVE 2008) these reactions are also strongly influenced by thermal 
conditions. 


In the present contribution nymphal growth of insects is presented in terms of a multi- 
compartment model including numerous aspects described above. The transfer rates 
between the compartments represent non-linear functions which among other depend on 
several environmental factors. The approach is used for the simulation of nymphal 
growth of the Australian field cricket Teleogryllus commodus. Thereby, theoretical 
results are compared with experimental data. 


Material and Methods 


Description of the theoretical model 


Larval or nymphal growth of insects represents a multi-phase process with each phase 
constituting a developmental stage. For that reason a multi-compartment model was 
developed, where each compartment corresponds with a different instar. Hence, growth 
of each larval or nymphal stage can be investigated separately (Fig. 2). From a 
mathematical point of view, the theoretical approach is founded upon the basic equation 
dK"/dt = (Wi? + Wn": j") + KP, (1) 
with K, Wi, and Wh, respectively, denoting the body length of the developing insect, the 
average intrinsic growth rate, and the average hormone-associated growth rate. The 
factor j determines the extent of Wn and thus ranges from 0 to 1. The factor n indicates a 
specific larval or nymphal stage. Mathematical solution of equation (1) results in the 
expression 
K*(t) = Ko": exp[(Wi? + Wa" - j") + t], (2) 
where Ko" describes the initial body length of the n™ larval or nymphal stage. Average 
growth occurring during a specific developmental stage can be expressed by the rather 
simple formula 
Wit + Wi -j° = [In(Ke") - In(Ko")]/L*, (3) 
with Kr" and L”, respectively, denoting the final body length and the temporal duration 
of the n™ larval or nymphal stage. The factor j" indicating the share of of hormonal 
growth in total size increase of the insect follows the equation 
J? = Jmin® - ((D" - D" - D®- L*) /(L”- L"/4)] - z, (4) 
where jmin denotes the minimal value of j, whereas the factor z constitutes the range of j. 


The variable D represents a specific time point within the phase of larval or nymphal 
development. 


858 


Total insect body length (K*") may be evaluated by using the final expression 
K' = Ko! - exp[(Wi! + Wai! - j!)- L'] + {Ko - exp[(W? + Wi? - j?) - L?] - Ko?} + ..... 
+ {Ko™- exp[(Wi™ + Wa” j™) > L™] - Ko™}. (5) 
In equation (5), L' - L™ denote the time spans of single larval or nymphal stages, whilst 
m represents the number of instars occurring during larvogenesis. 
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Fig. 2: Scheme depicting the multi-compartment model with its different parts. Number of 
compartments depends on the amount of developmental stages during larvogenesis. 


Experimental studies 


Nymphal growth of Teleogryllus commodus was investigated for environmental 
temperatures of 23°C, 27°C, and 30°C, respectively. Rearing of the insects was carried 
out under standard conditions (photoperiod: 12 h, relative humidity: 60%) in a speci- 
fically equipped climate chamber at the former Institute of Zoology, University of 
Salzburg (STURM 1999, 2002b, 2002c, 2003a, 2003b, 2008a, 2008b, 2011, 2012, 2014, 
2016a, 2016b, 2017, STURM & POHLHAMMER 2000). Small nymphal stages (3™ and 4" 
instar) were kept in small plastic boxes covered with a fine net, whereas larger stages (5" 
instar and older) were transferred into larger plastic boxes (50 x 30 x 20 cm) filled with a 
3 cm thick layer of peat soil. The boxes were equipped with egg cartons providing 
shelter for the animals. Nutriment of the crickets included fresh lettuce, porridge oats, 
standard diet for laboratory animals (Altromin 1222), and water. 


For an appropriate measurement and statistical evaluation of nymphal body length each 
larval generation was marked with a color code. Daily analysis of body length was 
carried out by selecting 5 individuals of each generation. The animals were anaesthetized 
in a CO: stream and measured from the forehead to the very end of the abdomen by 
using a slide caliper (accuracy: 0.1 mm). Antennae as well as appendices of the abdomen 
were intentionally excluded from the measurement. 
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Fig. 3: Comparison of theoretical and experimental temperature-development data of the black 
field cricket Teleogryllus commodus: (a) 23°C, (b) 27°C, (c) 30°C. 
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Results 


Time-dependent growth data obtained from theoretical simulations and related 
experiments are summarized in Fig. 3. In the laboratory studies nymphal development of 
Teleogryllus commodus was investigated for constant environmental temperatures of 
23°C, 27°C, and 30°C, respectively, whereas other external factors were kept constant 
(photoperiod: 12 h, relative humidity: 60%). Independent of the experimental line 
animals were supplied with standard nutriment (see previous section). In addition, a 
constant number of crickets was reared in the plastic boxes in order to minimize any side 
effects raised by intraspecific competition. 


As can be recognized from the single graphs depicted in Fig. 3, experiments conclude a 
reduction of the duration of nymphal development by more than 40%, when environmen- 
tal temperature is increased from 23°C to 30°C. Based on the high number of data points 
plotted in the diagrams and subsequent regressive analyses average nymphal growth rate 
can be quantified with 0.2 mm/d at 23°C but with 0.39 mm/d at 30°C. This corresponds 
to an enhancement of the growth rate by nearly 100%. Nymphal growth behaviour 
predicted by the multi-compartment model is characterized by slightly sigmoidal 
characteristic of the respective functions. This phenomenon, however, indicates slower 
growth of the insects at the beginning and the end of larvogenesis but faster growth at 
intermediate stages. Correspondence between theoretical predictions and results obtained 
from the laboratory may be evaluated as partly excellent. Depending on the 
environmental temperature, numerical discrepancies between model and experiment are 
on the order of a few percent. This circumstance results in Pearson’s correlation 
coefficients (r) amounting to 0.98 (23°C), 0.99 (27°C), and 0.99 (30°C). 


Discussion and Conclusions 


Larval and nymphal development of holometabolous and hemimetabolous insects is 
commonly characterized by the circumstance that a differentiation between intrinsic and 
hormonal growth can be made (BRYANT 1984, NIJHOUT & AMLEN 1998, CONLON & 
RAFF 1999, STERN 2003). Among insects with incomplete development nymphal growth 
never abandons completely as a consequence of deficiency or deprivation of nutrients, 
but its velocity may undergo a considerable cut-back (STERN 2003, STURM 2016a, 
2016b). Previous studies could clearly demonstrate that hormonal growth control nay be 
subject to temporary fluctuations (DAMOS & SAVOPOULOU-SOULTANI 2008, 2012, 
STURM, 2016), so that application of a variable parameter j introduced above seems to be 
highly plausible. Variation of this parameter with time allows to trace nymphal develop- 
ment of Orthoptera in the best way (STURM 201 6a, 2016b). 


In the case of the Australian field cricket correspondence between theoretical and 
experimental growth data may be evaluated as excellent for the most part. This impres- 
sion is further underlined by the circumstance that theoretical approach and experiment 
are responsive to eventual changes of the environmental temperature in a very similar 
fashion. For both research accesses any enhancement of temperature is associated with a 
gain of nymphal growth velocity. In contrast to previous approaches outlined by STURM 
(2002a, 2003a), the current model also considers a deceleration of nymphal growth after 
exceeding of an upper threshold temperature, thereby underlining its non-linearity. This 
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temperature exerting a negative effect on physiological processes is commonly situated 
between 34°C and 38°C (BEHRENS et al. 1983, HOFFMANN 1985, GEWECKE 1995, 
STURM 2002a, 2016a). 


The theoretical model presented in this contribution may be categorized as a mixture of 
empirical and biophysical approach, thereby largely considering non-linearity between 
environmental control factors and nymphal development. The biophysical aspect was 
realized by describing nymphal growth by means of intrinsic and hormonal growth rates. 
At this stage of model formulation temperature dependence of these growth rates was 
expressed by application of simple exponential equations (STURM 2016a, 2016b), where- 
as related concepts of VAN’T HOFF and ARRHENIUS were not implemented yet. A main 
reason for this simplification is represented by the circumstance that experimental data 
describing exactly the enzymatic processes in cricket nymphs are characterized by very 
limited availability hitherto. In addition, hormonal processes controlling the development 
within specific nymphal stages are not fully understood. Finally, establishment of a 
complex biophysical approach is not mandatory for answering specific questions of 
insect growth in a sufficient manner. 


From the content presented in this work it may be concluded that the nymphal growth 
model partly provides excellent results on the one side, but on the other side has to be 
understood as initial step with regard to the theoretical description of Orthopteran 
development. 


Zusammenfassung 


Theoretische und experimentelle Studie zum nymphalen Wachstum der australischen Feldgrille 
(Insecta: Orthoptera). — Der vorliegende Beitrag setzt sich mit einem Multikompartimentmodell 
auseinander, welches zur Simulation des nymphalen Wachstums von hemimetabolen Insekten 
befähigt. Generell lässt sich eine hormonelle Wachstumsrate (Wn) von einer intrinsischen (Wi) 
unterscheiden, wobei der zweite Faktor die Zunahme der nymphalen Größe aufgrund der 
Aufnahme von Nahrung inkludiert. Neben dieser grundsätzlichen Annahme berücksichtigt die 
theoretische Näherung auch den möglichen Einfluss von externen Faktoren (Umgebungs- 
temperatur, Nahrungsangebot, Tag-Nacht-Zyklus, intraspezifische Konkurrenz) auf das Insekten- 
wachstum. Anwendung und Validierung des Modells fanden für Teleogryllus commodus WALKER 
1869 statt, da bei dieser Spezies geeignete experimentelle Daten zur Nymphenentwicklung ver- 
fügbar sind. Vorläufige Modellergebnisse stimmen sehr gut mit vergleichbaren experimentellen 
Resultaten überein (r: 0,97-0,99). Grundsätzlich kann eine signifikante Abhängigkeit des nympha- 
len Wachstums von der Umgebungstemperatur festgestellt werden. 
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